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Abstract 

We study the spin 1/2 and spin 3/2 fermion fields in a thick braneworld scenario in six dimensions called 
string-cigar model. This smooth string-like model has a source that satisfies the dominant energy condition 
and undergoes a Ricci flow. We propose a new coupling for the fermions with a background gauge field 
which allows a smooth and normalized massless mode in the brane with positive tension. By numerical 
methods the mass spectrum and the massive eigenfunctions are obtained. The Kaluza-Klein massive tower 
exhibits the usual increasing pattern and, in this scenario, the coupling term does not allow tachyonic 
Kaluza-Klein states. The brane core and the background gauge field alter the properties of the massive KK 
tower, enhancing the amplitude of the massive states near the origin and changing the properties of the 
analogue Schrodinger potential. Furthermore, we find massive modes as resonant states in this scenario for 
both fermionic fields. 
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I. INTRODUCTION 


In the last years, the extra dimensions physics has acquired a conspicuous prominency, mainly 
due to the seminal models of large extra dimensions (ADD model) [T] and the warped compactifica- 
tion models (RS model) [2, j3|. The RS model, for instance, brought to light again the possibility of 
noncompact extra dimension [3], These braneworld models provide solutions for important prob¬ 
lems of the High Energy Physics, as the gauge hierarchy problem ra, the cosmological constant 
problem [5] and the origin of dark matter |6]. 

Fields endowed with a bulk dynamics provide an explanation for the brane stability by assuming 
a stable source configuration, as a topological defect [7j. In 5D (or more generically in codimension 
1) models, a domain wall brane is generated by a spontaneous symmetry breaking mechanism 
which also provides a width and inner dynamics to the brane [7]. Further, for bulk propagating 
fields reproduce the 4D physics, a massless Kaluza-Klein (KK) mode ought to have a compact 
support (localized) around the brane, what leads to a massless effective 4D action 1 9] . Nonetheless, 
unlike the gravitational and scalar fields, the vector gauge field has no localized massless mode 
j9], being required a coupling with the dilaton field [9]. For the spin 1/2 and 3/2 fermions, the 
localization of the massless mode also require an additional coupling, usually a Yukawa coupling 
EMU- The bulk dynamics and the additional coupling also provides the existence of KK massive 
resonant modes which brings important phenomenological consequences pME]. Another approach 
to obtain a normalizable massless mode for the gauge and fermionic fields in brane models with 
positive tension is accomplished by changing of the geometry structure for a Weyl geometry [24j . 

In six dimensions (or codimension 2), stationary braneworld scenarios with axial symmetry are 
called string-like models [261 [28] . The string-like branes inherited this name due their resemblance 
with topological defects in (3 + 1), as the cosmic strings |25|. Assuming a warped geometry with 
regularity conditions, the gravity is localized in the string-like models whose source is a global defect 
|251, an infinitely thin brane [28] or a local vortex [30]. The gauge vector field is also localized 
without any other coupling but the minimal gravitational coupling (311 - 133) . 

An important string-like scenario is the Gherghetta-Shaposhnikov (GS) model, whose bulk ge¬ 
ometry is the warped product between a thin brane and the two dimensional disk [28] . This vacuum 
solution of the Einstein equation with a negative cosmological constant traps gravity and provides 
a smaller correction to the Newtonian potential compared to the RS model |28]. Nonetheless, 
the GS model does not satisfy the regularities conditions at the origin and the dominant energy 
condition, as well m- 
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Giovannini et al found numerically a string-like solution wherein the source is a local Abelian 
vortex in the Einstein-Maxwell-Higgs model [30] . This smooth and thick string-like brane satisfies 
all the energy conditions and all the regularity conditions m- The components of the stress- 
energy tensor are concentrated around the origin for the winding number n = 0. For higher 
winding numbers, the source is shifted from the origin m- Yet, only a numerical solution for this 
model is known. 

String-like models with gravitational higher derivatives [41] were proposed and their cosmo¬ 
logical features analyzed [42] . Braneworld models in 6D was also proposed assuming non-trivial 
transverse manifolds. Carlos and Moreno found a smooth string-like solution with a cigar-like shape 
[35 ] whereas Kehagias proposed a conical tear drop-like transverse space to solve the cosmological 
constant problem [44] , Silva and Almeida used a section of the resolved conifold to build a thick 
string-like brane that traps the gravitational [45] . scalar and gauge fields |46]|. Other transverse 
manifolds include general Einstein spaces m, apple-shape [48], football shape [49], and others 

m- 

An analytical interior and exterior string-like model was proposed wherein the transverse man¬ 
ifold has a cigar-like shape, the so-called string-cigar model m ■ The bulk geometry is a warped 
product between a 3-brane and the so-called cigar soliton, a stationary solution of the Ricci which 
asymptotically converges to the 2D disk [52] . The Ricci flow appears in string theory as a RG flow 
of the worldsheet [53] and this geometrical flux has also applications topological massive gravity 
[ 54 ] and condensed matter m- 

The cigar shape provides a metric similar to that found numerically by Giovannini et al 130], 
which makes the metric satisfy all the regularity conditions at the origin and behave similar to 
the GS model at large distances m- Then, the string-cigar model can be regarded as a smooth 
extension of the GS model m- Since a Ricci soliton is an extension of the Einstein spaces, the 
string-cigar model generalizes the 6D model in Ref. m■ The Ricci flow in the transverse space 
makes the source undergo a flux that varies the brane-tensions and the bulk cosmological constant 
m- As a result, the source satisfies all the energy conditions and it has a bell-shaped core displaced 
from the origin, as in the Ref. I3D|. 

The spin 1 and spin 2 bulk fields have been already studied in the string-cigar brackground 
j56j 157 ]. The massless modes are localized at the brane core and they recover the usual string-like 
behavior at large distances EH El]. The Kaluza-Klein spectra are attained, showing an increasing 
behavior for both fields, and the correspondent massive modes are enhanced near the brane core 
mm- Besides, resonant massive states are present in the gravitational case [ 561 . 
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A relevant question to address is how to include the matter fields (fermions) in the thin and 
thick string-like models. Besides the spin 1/2 spinor which describes the ordinary matter, another 
important fermion is the spin 3/2 gravitino field (superpartner of the graviton) that arises in 
supergravity context and is a dark matter candidate [58]. Oda studied the localization of various 
spin fields in a thin string-like scenario |3U [32]. As the vector gauge field in 5D, the free spin 1/2 
and 3/2 fermions can only be localized (without additional interactions) on a string-like defect with 
the exponentially increasing warp factor (negative tension) [32]. Xiao Liu et al proposed to couple 
the fermions with a U( 1) background gauge field by which the zero mode is confined in the thin 
string-like brane with positive tension [33]. Parameswaran et al analyzed the massless and massive 
spectrum in a 6D supergravity model which enables a finite mass gap even for an infinite extra 
dimension [59]. Dantas et al [60] obtained a normalized massless mode for a fermion in a string-like 
brane with a transverse resolved conifold [45] . The depth of the potential well and the high of the 
potential barrier evolves with the resolution parameter [60]. However, although normalizable, the 
massless mode is not well-defined at the origin. 

In this work, we propose a new coupling for the fermions with a background gauge field which 
localizes the zero mode in the string-cigar model. The thin-string limit is considered, as well. 
Imposing suitable boundary conditions to guarantee the self-adjointness of the spinor operators, 
a normalized and everywhere well-defined massless mode is obtained for both the thin string and 
string-cigar models. As for the gravitational [56] and the U( 1) [57J vector fields, we find that the 
spin 1/2 field massless mode is shifted from the origin and sets around the displaced brane core. 
The mass spectrum is the same for both right-handed and left-handed chiralities and it is free of 
tachyons. Further, the KK spectrum has an increasing pattern which exhibits the usual linear 
behaviour of the Kaluza-Klein theories. For the Rarita-Schwinger field (spin 3/2), the zero mode 
and massive spectrum have minor changes when compared to those of the spin 1/2. However, the 
amplitudes of the massive eigenfunctions for the spin 3/2 are higher than those for the spin 1/2. In 
comparison with the thin string-like model, the core of the string-cigar brane enhances the massive 
modes near the origin for both spin 1/2 and 3/2 fields. Besides, in a Schrodinger approach, the 
spin 1/2 and 3/2 fields possess identical behaviour. Furthermore, the coupling also enables the 
presence of resonant modes (massive modes solutions having very large amplitude near the brane). 

This paper is organized as follows: in section [IT] we review the features of the thin string and 
the string-cigar models. We comment on the main properties of the gravity and gauge fields in 
these models, which were developed in References [56] and m , respectively. We also studied the 
dynamics of the scalar field in these scenarios and show that is it identical to the gravitational case 
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56]. In sections III and IV, we study the localization of spin 1/2 and 3/2 fields. Furthermore, 


the massive spectrum is obtained in both cases. Concerning the resonant modes, we conclude 
that the Schrodinger-like potentials are the same for both spin 1/2 and spin 3/2 cases. Thus the 
behaviour of the resonant modes are only studied in subsection |IIID| In section [Vj conclusions and 
perspectives are outlined. 


II. THE STRING-CIGAR BRANEWORLD 

Consider a six dimensional spacetime Mg where the 3-brane M 4 can be embedded. Assuming 
an axially symmetric and static bulk Mg, the braneworld scenario is called a string-like model 
|26l - t28l [30] . A general metric for the string-like model takes the form [2614281 (30lj 

dsg = F(r)r] fiu dx^dx 1 ' + dr 2 + H(r)d8 2 , (1) 

where 0 < r < r max and 6 e [0, 2n) are the radial and the angular coordinates, respectively. The 
radial component can extend to infinity, i.e., ?’ max —> 00 . In order to guarantee that the scalar 
and the extrinsic curvatures are finite at the origin, it is usual to impose the regularity conditions 


F(0) = (VW))' r=0 = 1, and F'(0) = H(0) = 0, (2) 

where the primes denote derivatives with respect to r. 

Performing the change of coordinate [29] 

z (r) = f F~^ (r')dr', (3) 

Jo 

the metric (jT]) can be cast in the conformal form [29] 

dsQ = F(z) (rj^dx^dx 11 + dz 2 + /5(z)dd 2 ) , (4) 

where /3(z) := fgj. 

Let us assume that the bulk dynamic is governed by the Einstein-Hilbert action with bulk 
cosmological constant A mm Enj: 

Sg = f (--^R - A + yJ^gdPx, (5) 

where kq = 8tt/Mq , Mg is the six-dimensional bulk Planck mass and L m is the matter Lagrangian 
for the source of the geometry. From the matter Lagrangian L m we define the stress-energy tensor 

T = T MN dx M ®dx N , (6) 
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where the stress-energy tensor components are defined by [26l - f28l [30] . 


Tmn = 


2 d(^gL r 


V=9 dg MN ' 

The Einstein-Hilbert action provides the Einstein equation for the bulk 

Rmn — ^ r 9mn = —K§(A.gMN + Tmn)- 

An axisymmetric and static anstaz for the stress-energy tensor has the form EMSGSJ]. 

T = to(r) ^e 0 0 e 0 + 0 + t r (r)e r 0 e r + 0 e^. 

Using the metric ansatz (jTJ) and the stress-energy ansatz the bulk Einstein equation ( 
to the system of equations 

3 f F'\' 3 (F’\ 2 3 F'H' 1 (H'\ 2 1 (H' 

2V^y + 2 V^fJ + aT~H + 1\~H 


2 V H 


3 /ry F_FV_ 
2\F) + F H 


F' 


!\ ' 


5 ( F' 


2 V F 


( 7 ) 


( 8 ) 


(9) 

yields 


K G (A + t 0 (r)), 

( 10 ) 

k 6 (A + t p (r)), 

( 11 ) 

k 6 (A + t e (r)), 

( 12 ) 


A. Thin string-like model 

In the Gherghetta-Shaposhnikov (GS) model, a vacuum solution of the system of Einstein 


equations ( 10 ) was found 


Assuming that -y = -y = 0, the metric functions are given by 


F(r) = e~ cr , H(r) = B%F(r), 


(13) 


where Rq is an arbitrary length scale and the constant c is related to the bulk cosmological constant 
by [28] 


c 2 = —7 k 6 A. 


(14) 


Hence, the cosmological constant must be negative and the bulk is an AcISq spacetime (28]. Since 


the functions in Eq. (13) are vacuum solutions, the GS model represents an infinitely thin string-like 


braneworld which is an extension of the Randall-Sundrum (RS) metric to six dimensions. 

For the thin string-like metric, i.e., F(r) = e~ cr and H(r) = RgE(r), the conformal coordinate 
can be found to be [29] 


z(r ) = -|e2 —l) . 


(15) 
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Note that z(r = 0) = 0 and z'(r) > 0 in Eq. (15). Then, the conformal coordinate z(r) is still a 


gaussian radial coordinate which measures the distance from a point in the transverse manifold to 


the origin. Using the conformal coordinate z in Eq. (15), the metric factors have the form 

4 1 


F(z) = 


2\ 2 


(*+l) 


H(z) = RqF(z). 


(16) 


Furthermore, the hierarchy problem between the four-dimensional Planck mass (M 4 ) and the 
bulk Planck mass (M 6 ) is solved in this scenario [28] [51], and these masses are related by the 
following equation 


Ml = 2vrAf 6 4 / ^-g{r)F~ l {r')dr' 

Jr'=0 

2lT Rq f s nix 

= —l 1 " 6 * 


(17) 


The relation between the bulk and brane Planck energies can be rewritten using only the ratio 
between the bulk cosmological constant and the string tension [28] . 

For fluctuations of the metric (jTj) in the form [28], 30, 5lj 

ds\ = F(r)(rjn„ + h^ w )dx+dx u + dr 2 + Ff(r)d9 2 , (18) 

satisfying the traceless gauge V / ''h /i „ = 0, the linearization of the Einstein equations ([8]) yields the 
equation for the gravitational perturbation [28] [51] 

AB ‘ 


AlO 


(19) 
[51] and 

imposing the free wave dependence on the 3—brane □^^(x) = m 2 /r Atl/ (x) [28, 00] El], the radial 
component of the graviton equation of motion takes the general form gg G 2 H EU 


= 9a(V- 96 V d B hfj, v ) = 0 . 

Performing the Kaluza-Klein decomposition /^(x, r, 6) = h lw {+) Y1 Ai=o 9n,l( r ) 


d 2 + 27 d r + 


f m l,l 

?x\ 

u 



<i>nA r ) = °> 


where 


m/ . F' 1 H' 5 F' 1 p 
^ r ~ F + 4H “ 4 F + 4 /T 

The graviton radial equation for the thin string-like model has the explicit form 

Am - \ C( Am + { m 0 ~ l 2 /Rl) e ° r = 0 . 

From the Eq. Q, we impose the boundary conditions gSEQlEI] 

4> (oo) = (f>' (oo) = 0 . 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


7 












The radial graviton equation (22) together with the boundary conditions (23) forms a Sturm- 


Liouville problem |281130] 5lJ. For m = 0 and s-wave solution {l = 0), the localized massless mode 
solution is obtained as 


4 > m= o(t) — 


3c 

2Rq 


— ~cr 

e 4 cr , 


(24) 


~ 3 

where <\> m = e“4 cr <fr m=0 


For m/0, the massive modes has the form 


4>m{p) = e-i cp 




+ B2Y5/2 I ez 


2m 1 


cp 


(25) 


where B\ and B 2 are arbitrary constants and m = m'^ — P/R^. The exponential dependence reveals 


that the massive modes are not localized on the brane [28]. Applying the boundary conditions (23) 


on the massive modes (25), the graviton mass spectrum in the GS model was found as 


m r 


vr _ £ 

n -— e 2 

2 2 


(26) 


where r max is a finite radial distance cutoff. The gravitational massless mode is localized in the 
brane and the contribution from the nonzero modes provides a small correction to the Newtons 
law on the 3—brane [28j . 

The vector gauge field was also studied in the string-like models [3IH331 136) [37] 46] [57] [59] . 
Starting with action 


Ss P in -1 = f V=gg MN g RS ?MN?R S d 6 x, ( 27 ) 

where ‘Jmn = V m-A-n — ^nAm, the equation of motion is [3TTf33] 


J_ g {d s V=~gg SM g RN ?MN) = 0 . (28) 


Imposing the gauge conditions [[3TH33] d^A^ = Ag = 0 and A r = A r (r,6), 

read [31] [32] [Ml EZ] 

the Maxwell equations 

[g^ u d^d u + A r = 0, 

(29) 

/ F 2 \ 


o' 

II 

(30) 

[rrdf/A + ^ d 2 e + -±=d r (FVHd r ) j A a = 0 . 

(31) 


( JS 

Using the Kaluza-Klein decompositions A^{x M , r,9) = M A.)?' \x fl )p n (r ) e lW and A r (x M ,r, 9) = 

n,l=0 

OO / 7\ 

E (x^)g(r) e lW , the radial dependence of the gauge field on the brane for l = 0 is governed by 

1=0 
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the Sturm-Liouville equation [46] [57 1 


F' 


%+[2v- — )a r + -£ 


mz. 


Pn(r ) = o, 


which can be rewritten as 


3 F' l/3'\ , 


mt 


Pn(r) + ( - T + ) p'n(r) + Pn(r) = 0, 


(32) 


(33) 


where f3(r, c ) = H(r, c)/F{r, c ). For the thin-string model, the radial equation (33) has the explicit 
form ED (323 


4>m - 2 + ( m 0 - l 2 /Ro) e ° P Pm = 0. 


(34) 


Unlike the vector gauge field in the RS model, in the thin string-like model the gauge field has a 
massless mode of form mm 


Po(r) = 


5c 
2 R~o 


e-F, 


which is normalizable ED|32]. The 
Pm(r) = e2 f 


massive 


modes have the form mm 


2m 


BiJ 3/ 2 — ei cr + B 2 Y 3/2 — eT 


2 m 


(35) 


(36) 


and then, are non-localizable ED [32]. Imposing the same boundary conditions ( |23| ), a linearly 
increasing KK gauge spectrum was found ED EH 


m n = —nire 2 


Finally, for a minimally coupled scalar held, the action ED EH 05] 

Sspin-o = -7, J \f^gg MN d M $d N <&d 6 x, 
provides the equation of motion 

1 




d M [V=lig MN d N ^>] = 0 . 


(37) 


(38) 


(39) 


Using the general string-like metric (jT]) in Eq. (39), we have 

( 


, n ,w dr(F 2 H 2 dr 

F M F 2 H \ 


+ « 

H 


$ = 0 . 


(40) 


Using the KK decomposition <b(x M ,r, 6) = -^^(a^) D Xn{r)e lie and r]^ u d fl d u ip(x tl ) = m 2 tp(x^), 

n,l =0 

the Eq. (39) yields the radial sturm-liouville equation 

(TTi^ 72 

d 2 r +2y{r)d r + l ) X n(r) = 0. (41) 


f m l,l 

l2 \] 

u 

»)\ 
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Note that this equation is the same for the gravitational case (20), thus for s-wave solution, both 


zero mode and as the KK spectrum are given by the results of the graviton (321051 EH EB!. This 
behavior between Spin 0 and spin 2 modes is also verified in five dimensional scenarios [5]- 

Despite the good results of the thin string-like model described above, the regularity conditions 
Q at the origin and the dominant energy conditions are not satisfied for this model |30l [40] . This 
issue arises due to the metric be only an exterior solution of the Einstein equations. Although it 
is possible to consider an interior and exterior solution separately using the junctions conditions 
to match the solutions, the thin string limit (where the width of the core vanishes and only the 
exterior solution remains) is known to present the same problems EDI- As a matter of fact, 
Giovannini et al have found numerically an interior and exterior string-like solution of the Einstein 
equations satisfying all the regularity and energy conditons whose source is an Abelian vortex in 
the Einstein-Maxwell-Higgs model. Nonetheless, the analytical expression of the solution is yet 
unknown. 


B. String-cigar model 

An analytical interior and exterior smooth extension of the GS model, called the string-cigar 
model, was proposed in Ref. EH • The thin string-like GS model is built from the warped product 
between the 3—brane and two dimensional disc O 2 of radius Rq, whose metric has the form 


dsi, 2 = dr 2 + RidO 1 


(42) 


The constancy of the radius prevent the GS model to satisfy the H( 0) = 0 and = 1 

conditions. 

In the Ref. EH, a regular transverse manifold C 2 with metric 


1 


dsi = dr 2 + — tanh 2 (cr)dd 


(43) 


was proposed. Note that asymptotically, the metric (43) converges to the D 2 metric (42), for 


Ro = 1/c. At the origin, the effective radius tanh(cr)/c vanishes. Hence, the transverse manifold 
C 2 has a cigar-like behaviour. 

Indeed, the transverse space C 2 is a stationary solution of the geometric Ricci flow according to 


dc/ab 

dc 


= -2 Rab{c), 


(44) 
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where c can be regarded as a metric parameter and R a b is the Ricci tensor. The metric solution 


(43) is a Ricci flow solution called cigar soliton [5BE2]. The Ricci flow ( |44[ ) defines a family of 
smooth geometries depending on the parameter c m- Using the cigar soliton 62 as the transverse 
manifold, a string-like model was proposed as m 

F(r) = e - [cr - tanh(cr) ], (45) 

and 

H(r) = (3(r)F(r), (46) 

where /3(r) = tanh 2 (cr). Since this string-like mode has a cigar transverse manifold, it was 
named string-cigar model. We plot in the Figure [l] the metric functions for the thin string-like and 
for the string-cigar models. It can be seen that, for large r, both models possess similar behaviours. 
Therefore, the string-cigar model recovers the thin string-like limit asymptotically. At the origin, 
the string-cigar model satisfies all the regularity conditions m- 


From the Einstein equations (10), the components of the stress-energy tensor corresponding to 


(47) 

(48) 

(49) 


the string-cigar braneworld are m 

c 2 / 13 5 \ 

f 0 (r, c) = — (7 sech 2 cr -1-sech 2 cr t anh cr -sech 4 cr) 

ACg \ 2 2 / 

c 2 / 5 \ 

t r (r, c) = — (5sech 2 cr + 2sech 2 cr tanh cr-sech 4 cr ) 

Atg V 2 / 

c 2 / 5 x 

f@(r, c) = —( 5sech 2 cr + 4sech 2 cr tanh cr-sech 4 cr). 

Atg V 2 / 

The components are non-negative to,t r ,tg > 0 and satisfy the dominant energy condition, i.e., 

to > |t r |, We plot in the Figure [2] the energy density to for different values of the parameter c. 

The position of the maximum of the energy-density evinces that the brane core is shifted from the 

origin, as for the Abelian vortex model for higher winding number m- 

Since the stress-energy tensor components have a compact support around the origin, the string- 

cigar geometry can be regarded as an interior and exterior thick string-like solution of the Einstein 


equations (10). Once the source evolves with the Ricci flow parameter c, the string-cigar model 
reflects changes that both the brane source and the bulk cosmological constant undergo [51]. For 
instance, due to the tanh(cr) and sech(cr) terms, the stress-energy components vanish when c —> 00 
and we recover the thin GS model. 

In the string-cigar model, the relation between the bulk and the brane energy scales, given by 
Eq. pT] ), reads (HJ 

2ttMr I ~ e 2' 


Ml = 


f 


00 — | (cr—tanh(cr)) 


t anh (cr) 


-dr. 


(50) 
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Then, in order to guarantee that M 4 Mq, the parameter c must be small m- 



r 


Figure 1: Warp factors for c = 0.9. The thick lines correspond to the string-cigar geometry, while the thin 
lines to the GS model with R 0 = 8.0. 



r 


Figure 2: Energy density of the String-Cigar braneworld for different values of c. Its maximum indicates 
the position of the core of the brane. 


The localization of gravity and U{ 1) vector gauge held in the string-cigar braneworld was 
performed in References 1SU and m , respectively. 

The radial graviton equation has the form [51] 


+ c 


5 , 9 . . sech 2 (cp) 

— tanh 2 (cp) + -—— 7- 

2 v tanh (cp) 


0' m + («>» (m 2 n - An = 0. (51) 


Note that asymptotically, tanh(cr) = 1 and sech(cr) = 0, then the graviton equation in the string- 
cigar model has the same form of the thin string-like GS model [51]. Nevertheless, near the origin, 
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the effects of the brane core change the behaviour of the gravitons. In fact, its massless mode is 
given by m 


4>m=o(r ) = NiFi(r)^a(r), 


(52) 


where N\ is a normalization constant. The massless graviton mode (52) satisfies the boundary 


conditions (23) and asymptotically behaves as the GS massless mode. Further, the zero mode has 
its peak shifted from the origin, as the energy density [51]. 


The massive modes as solutions of the Eq. (20) was studied in details in Ref. m , where the 
complete Kaluza-Klein spectrum and the corresponding eigenfunctions were attained. Asymptot¬ 
ically, the massive KK states has a similar behavior of the thin GS model whilst near the origin, 
the brane core interaction enhance their amplitude compared with in the GS model. The resonant 
states were also found as solutions of the analogue Schrodinger equation [56]. In the string-cigar 
model, the graviton massive spectrum, obtained numerically, has a linear behaviour, as for the GS 
model [55] . 

Since the radial equation for the scalar field is the same of the graviton, the analysis of the 
massive modes for the scalar field in the string-cigar braneworld is identical to the gravitational 
field presented in Ref. [SS3- 

For the vector gauge held, the radial equation reads m 


Pn T c 


3 2 sech 2 cr 

— tanh cr -\ --- 

2 tanh cr 


^ + e (—tanhc,) m 2 pn = 0 . 


(53) 


Asymptotically, the gauge KK equation (53) has the same form of that for the thin string-like 


model. The localized massless solution of the Eq. (53) in the string-cigar scenario were found in 
Ref. [57] as 

Pm=o(r ) = N 2 F^ (r)P J (r), (54) 

where N 2 is a normalization constant. It is interesting to compare the massless mode of the bosonic 


and fermionic fields. We insert a comparative plot of massless modes of gravity, Eq.(52), vector 


held , (54) and fermionic helds in the Sec. IV 


The massive solutions of Eq. (32) for the string-cigar was studied in Ref. m- The well-known 
linear increasing behavior for m <C c was obtained. As for the gravitational case [56], the massive 
eigenfunctions behave as the thin string-like modes asymptotically and they are influenced by the 
core of the brane near the origin. However, resonant states were not found m- 

The similarity between Kaluza-Klein modes of the spin 0 and spin 2 helds is expected. This 
behaviour is present in the thin string-like model as inferred by Oda in Ref. [32]. Besides, the 
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same feature occurs for the of spin 1/2 and spin 3/2 cases, which will be discussed in the sections 

Iffll and [El 


III. SPIN 1/2 FERMIONS 


Consider the action on curved space background for bulk massless spin ^ fermions 

5 6l/2 = [ V^mT M D M *<Fx, 


(55) 


where T M = are the curved Dirac matrices defined from the flat Dirac matrices T jW through 


the vielbeins gMN = CaHnVmn- These matrices obey the Clifford algebra {r JW ,r JV } = + 2 g MN ts. 


M rAb 


,ATJV 1 


Here, Dm is the gauge covariant derivative given by 


Dm = 9m + Dm — iq^-M, 


(56) 


where Dm = \v PN £n ^M^p + ^mq^p is the spin connection and Am = A fl (x)x + A e (r)8 

is a cylindrically symmetric background gauge vector field. We remark that here, Am is not related 
to the dynamical field Am of the previous section. 

The non-vanishing terms of the spin connection are 


1 F'ir) „ „ , n 1 H'(r) 

D u = - . = TnTf and Dq = - - 

M 4 jF[v) 11 4^ 


TgTj.. 


(57) 


Substituting the equations (56) and (57) in the action (55) in the background metric given by Eq. 


([!]), we obtain the Dirac equation: 


T m D m ^= 


F-ir^d, - iqApixj) + r(a r + y + ^)+ 


(58) 


T = 0. 


+ H— 2 T 9 (d e - iqA e {r) 

We will now choose the usual Weyl spinor and gamma matrices representation on six dimensional 
models [3211311351 EMU: 


\k(x, r, 6 ) = 


0 


(59) 




r = 


(60) 


0 7 | r = [ 0 75 

0 y 0 / \r 

Due to the metric sign convention (—, +, +, +, +, +) the 7 ^ matrices in Weyl representation become: 


0 -Y 
,e o 


7 0 = — ia 1 


1 ; 


7* = —a 2 <g> a\ 


7 5 = cr 3 (g> t 2 , 7 e = H 4 , 


( 61 ) 
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where u l are Pauli matrices and 1 is the identity matrix. In this convention, the matrix 7 0 is 
anti-hermitian, while the others are hermitian. The y 5 is such that 7 5 i^ R ,l = The Dirac 

operator acts as 7 ^ ( d^ — iqA Others representations can be directly deduced from the 
general forms of the Ref. [62] for signature (+, —, —, —, —, —). 

Now, let us perform a Kaluza-Klein decomposition on t/q in the form 


Mx, r, 9) = ~^= Y [^R nA {x)a Rn l (r) + ^ Ln l {x)a Ln l (r) 

v n,l 


Aid 


(62) 


Using the equations (59), (60) and (62) for the s-wave solution mmm, the Dirac equation 


(58) turns to the following chiral coupled equations 


[d r + T(r) + W(r)] a Rn (r) = --^js =a Ln (r) 
[d r + (P(r) - W(r)] a Ln (r) = -^p.a Rn (r), 

where 


(63) 


and 


IP(r) 


F'(r) ( H'{r) 
F(r) + 4 H{r) 


—c 


— tanh 2 (cr) — 2 sech 2 (cr) coth (cr) 


(64) 


W(r) = q 


A e (r) 

VW) 


cq 


Ae(r) 
tanh (cr) 


e |[c 7 —tanh (cr)] 


(65) 


A. Spin 1/2 Zero Mode 


For to = 0, the expressions (63) decouple in two first order differential equations which solutions 


are: 


a Rn,L s (. r ) = Co exp 


- / (3>±W )dr' 


( 66 ) 


where Co is a normalization constant. 

In order to the zero-mode be localized and free of singularities at the origin, we impose the 
orthogonality condition 


W Rn ,L s {r)\ dr' = 5 RntLa . 


(67) 


This condition implies that lim,.-^ L (r) = 0. Nonetheless, since 

5 


— / y{r)dr' = 


2 (tanh (cr) 

cr — tanh (cr) -|— in - 

5 V c 


( 68 ) 
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is non-convergent, the Ag(r) function presented in Eq. (65) has to be adjusted in order to overcome 
this drawback. Assuming that 


W (r) = -AT(r), 

where A is a dimensionless coupling constant, the rigth-handed zero mode becomes 


a R.(, r ) = Co exp 


dr' (A - 1) T 


= C 0 F (A - 1) (r)F3 (A - 1) (r). 


(69) 


(70) 


For A = 0 this expression is the same obtained in Ref. [32] which is non-normalizable. 


Using the explicit expressions of the warp factors in Eq. (45) for string-cigar, the zero mode 
and the gauge angular component are given, respectively, by 

/ tanh (cr) 


2 ex p ( 


(A-l) 

2 / 5 

exp ( -(1 — A)[cr + tanh (cr)] 


and 


A e (r) = - 
Q 


5 

- tanh 3 (cr) — 2 secli 2 (cr) 


| [cr —tanh (cr)] 


(71) 


(72) 


In the absence of the coupling (A = 0), the massless mode is not localized in the brane. For A > 1, 
the zero mode is normalizable, but only for A > 3 its derivative is continuous and null at the origin. 
From these restrictions over A, we find the following boundary conditions 


a 


d r 


%l(P) = Uma% L {r) = 0 


(73) 


a °R,L(r)\ r=o = r lhn d r [a% L (r)] = 0. 

We plot the fermionic zero mode and the gauge angular term in the figures [3] and [4] respectively. 
In both cases, the coupling constant A controls the amplitude, whereas the geometric parameter 
c regulates their distribution over the radial extra-dimension. The displacement of the zero mode 
from the origin is an important result present in expression This feature is closely related to 

the fact that the brane core is not placed at r = o m ■ This zero mode has a similar profile to 
that of the energy density of the string-cigar model I5U- Note, in the figure [3j that the zero mode 


satisfies the homogeneous boundary conditions (73) provided that A > 3. 

In order to confine the massless left-handed fermions we would have to make A —> —A in the Eq. 


(69) with the same restriction |A|> 3. Therefore, only one massless chiral mode can be trapped in 
the brane. This is a well-known result in five warped dimensional models [12 [IIES [EMI]- 
It is interesting to note at this point that the references [32l, 33] use a less restrictive imposition 
on the radial component oP R L ( r ) in the form 


h = 

2 


poo poo 

/ Ii (r)dr' = / ^F- 2 (r)\a° RntLs (r)\ 2 dr' = 5 Rs , Ln , 

Jo 2 Jo 


(74) 
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Figure 3: Plot of the right-handed fermionic zero modeFigure 4: Plot of the gauge angular component for the 
for c = 0.5 in the String-Cigar model String-Cigar model. 


which comes from the effective action S^jj(x, r, 6) 

/ oo roc 

^(x)^ [d,j, + iqA^x)] ip(x)d 4 x' / Ii (r)dr' 

-oo Jo 2 


using the equations (55) and (62), namely 
d6 

o * Jo ^ir 


(75) 


However, if we adopt only the condition above, the spinor in the string-cirgar scenario will be 


exposed to singularities, while for the condition (67), these singularities vanish. Besides, the Eq. 


(74) is satisfied too. 


B. Spin 1/2 Massive Modes 


In order to study the massive modes, let us decouple the system of first-order differential equa¬ 
tions (63) performing the conformal change of variable z(r) given by Eq. © which turns the 


coupled first-order differential equations system (63|) to 
(d 2 + 2?8 Z + ? 2 - W 2 + (IP ± W 

where 


aR n ,L n (z) = -m 2 n a Rn , Lri (z), 


y = y(z)y/F{z) , w = w {z)J¥{z). 


(76) 


(77) 


Here, the over-dots mean derivatives with respect to z and the +, — sign stands for the right and 
the left chirality, respectively. 


The decoupled system of two second-order equations (76) is composed by independent Sturm- 


Liouville problems for each chirality. Thus, we can analyse the dynamics for the chiralities in¬ 


dependently. However, due to the involved form of the metric components (45), the conformal 
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coordinate z(r) in Eq. ([3]) can not be achieved analytically in general for the string-cigar geometry. 
The functions !P(z) and W(z) (and their derivatives) must be constructed from a numerical integral 
of Eq. ([3]). In order to avoid the cumulative round-off errors in the forthcoming analysis, we will 


study the Eq. (76) in the r coordinate, where the metric functions are already defined. It turns 


out that returning to the r coordinate, the second-order system (76) is still decoupled, and it can 
be written as 


a R n ,L n ( r )+ 3 / + 2 g a 'Rn,L n ( r ) + I ^ ^ ^ 


+ (1-A 2 ) 


5/ 2 + fg + 10/ + 2 sf 


+ 


i 2 


4 /+ 4J 


ml 


>a R n ,Ln = -ya B n , Ln (r) 


where 


\ F'(r) , , , P'(r) 

fir) = --and g(r) = 


Fir ) 




(78) 


(79) 


1. Thin string 


For the thin string-like model (/ = — c, and g = 0), the Sturm-Liouville KK equation (78) 
reduces to 

a R n ,L n i r )~ 3ca 'R,Li r ) + Tr[( 1: h A ) + _ x2 )] a Rn,Lnir) = ~m n e cr a^Lnir). (80) 


For A = 0 (absence of the coupling), the Sturm-Liouville equation (80) turns to 

35 

a R,Li r ) ~ 3c «ki( r ) + ^ c2 " R T( r ) = -m n e cr a RntLn ir), 
and the solutions have the form 
a R,L = Ar l e 


(81) 


3 cr 
2 


( 2m cr \ i f 2m cr 

J ±\ [~ e2 ) +B R,l y ±\ { — e2 


(82) 


where L , L are integration constants. It is worthwhile to mention that the massive modes 
depend on the Bessel functions of order nn,L = while the graviton has order 5/2 (25) [28] and 


the gauge vector field has order 3/2 (33) 


Moreover, unlike the graviton, gauge vector field and scalar field, the fermionic massless mode 


for A = 0 in not localized on the thin string-like brane (32]]. In fact, for m = 0 in (81), the massless 
solution has the form 


a R,L ~ A° r l ePlCr +Br,l e 


p 2 cr 


(83) 
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where p \.2 = and A Q R L , B R L are integration constants. 

For A 7 ^ 0, we write the solution as 


olr,l = Ar,l e 2 


•W • +B ^ 


( 2m cr 

^ (v e2 


where Arj^, Br^ are integration constants and 


Pr,l = (5A ± l)/2 


(84) 


(85) 


are the orders of the Bessel functions. The massive modes in Eq. (84) bears a resemblance with 
those found in 5D for massive fermions 


j, where the mass, as the gauge coupling A here, controls 
the order of the Bessel functions. We plot in the figure [5] the analytical right handed solution (84) 
for different values of A. Note that the gauge coupling distances the massive modes from the brane. 



Figure 5: Massive mode in the thin string-like (GS) scenario for m = 0.50. The values of the parameters 
was An = Bn = 1.5 x 10 -3 , A = 5.0 and c = 0.5. 


Unlike the massless mode, the massive eigenfunctions (84) are not trapped in the thin string¬ 
like brane due to the exponential and the Bessel functions. Nevertheless, in order to satisfy the 
boundary conditions, the order of the Bessel function ought to be pr > 7 and pl > 8 . The coupling 
allows the Bessel functions order to be integer or half-integer. For A even the order pr^l is half 
integer whereas for A odd pr,l is an integer. Although the coupling constant A can be any real 
number, the Bessel functions of irrational order suffers of branch issues and then, henceforward, 
we shall be concerned with the rational A only. A noteworthy feature is that the massive modes 
are related by pr = pl + 1. An interesting reason for this symmetry will be shown in the next 
section through the Schrodinger approach. 
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Applying the boundary conditions (73) in the massive modes (84) at the origin and at some 


cutoff distance r = r max , for m < c we obtain the conditions Brl = 0 and 


. 2 Tfl cr max . 

e 2 1=0- 


( 86 ) 


From the roots of the Bessel function (86), we find that the KK massive spectrum m n is discrete 
and it behaves as the series m 


m r 


C7T _ crmax 

~2 6 


. 2 fi R L - 3 hr l (2 - Hr,l) 

n H---1-— 


n 




+ o ( —2 

n z 


(87) 


The KK spectrum (87) exhibits an increasing behaviour, as expected ||32j . For large n, the spectrum 


behaves linearly whereas for small n the 0 (^) terms in the series (87) changes the rate of increasing 
of the masses. The mass gap between the massless and the first massive mode is given by 


. C7T _ c ^max 

A = m 0 « — e 2 


( 88 ) 


which vanishes for an infinite radial coordinate. Then, for an infinite radial extra dimensions, 
there is no mass gap between the massless mode and the massive KK tower, as usual in warped 
compactihed models 


2. String-cigar model 


For the string-cigar model, where 


f(r ) = —ctanh 2 (cr) and g(r) = 2c 


sech (cr) 
tanh (cr) ’ 


(89) 


the eigenvalue problem (78) is quite complex to be studied analytically. Then, we employ numerical 


methods to find the mass spectrum m n and the correspondent eigenfunctions. The numerical 


integration of the Eq. (78) was performed using the matrix method [63] based on finite differences 


with second order truncation error. To avoid the singularity in r = 0 and the overflow errors 
provided by the exponential functions for large r, we discretized the domain r e [0.01,13.00] into 
an uniform grid with constant stepsize h = 0.01. 

We plot the lowest mass eigenvalues in the figure [6] for A = 7.0 and different values of the 
geometric parameter c, which is related to the bulk Planck mass m- Note that the spectrum 
is monotonically increasing, as expected from the Kaluza-Klein theories. This regime is valid for 
m <C c m ■ Further, heavier masses will be acceptable as c increases. Moreover, the growth rate 


of m n is slightly lower for the first indexes n. This is in accordance with Eq. (87). It is worthwhile 
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to mention that, although we treat both chiralities left and right in Eq. (78) in an independent 


way, the relation among the eingenstates reveals that, regardless the massless mode, the left and 
right spectrum are the same. 



n 


Figure 6: Mass spectrum for A = 7.0 in the String-Cigar model. Note that heavier masses will be acceptable 
when the bulk Planck mass increases. Moreover, the growth rate is slightly lower for the first indexes. 


In the Figures [7] and [8j we present the eigenfunctions for c = 0.5 and for A = 5.0 and 9.0, 
respectively. Near the brane, they behave as Bessel functions of integer (> 2) order which increases 
with A. Since the string-cigar model recovers the thin-string one asymptotically, it is expected 
that the eigenfunctions have the same behaviour for large r |51j . Note that this occurs when 
compared with Figure [5] Moreover, the core of the string-cigar brane amplify the massive modes 
near the origin. This behaviour occurs for the gravitational m and gauge m fields. This is the 
first stimulus for searching resonant modes. In the next section, we will present the formalism 
concerning resonant states. 


C. Schrodinger approach 


The equality of right and left spectra and the relation between the eingenfunctions result from 
an underlying symmetry which is manifested in the Schrodinger approach. Indeed, performing the 
change of variable 


ur,l{z) = exp 


- / ?{z)dz' 




in the Eq. (76), we obtain a Schrodinger-like equation as 


(90) 


[-<9 z + Vr,l(z)\ olr )L (z) = m 2 a RtL (z), 


(91) 
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Figure 7: Normalized eigenfunctions for c = 0.5 
and A = 5.0 in the String-Cigar model. The 
masses eigenvalues was obtained as m R = 0.4024 
and rriL = 0.4025. 



Figure 8: Normalized eigenfunctions for c = 0.5 
and A = 9.0 in the String-Cigar model. The 
masses eigenvalues was obtained as m R = 0.4593 
and rriL = 0.4586. 


where 


v R , L (z) = w 2 (z)±d z w. 


(92) 


For the thin string-like brane, the analogue potential has the form 

1 


Vr,l( z ) = v 


— =F 1 
2 


(93) 


J ( z +l) 2 ' 

For A > 0.4, there is a potential barrier at the origin for both chiralities. Since the potential 
vanishes asymptotically, there is no mass gap for an infinite radial coordinate, as we found using 


the Sturm-Liouville approach in section III B Defining the variable 


x := m ( z 4— 


the Schrodinger-like equation for the thin string-like model reads 


-«+ T 


2 


jx 2 f 1 


0!R,l(x) = 0, 


whose solution can be written as 


(94) 


(95) 


&r,l(x) = N RtL y/x J^ L {x) + A R:L Y^ rl (x) 


(96) 


The Eq. (96) is only another form of the massive KK state given by Eq. (84) 
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The structure of the potential (92) enable us to rewrite the Schrodinger-like equation as the 


system of equations 

H R a R {z) = m 2 a R (z), H L a L (z) = m 2 a L (z ), 
where the Hamiltonians operator H R l can be factorized into 

H r = A* A, H l = AA+ 

and 

A(z) :=^+W (z). 

For the thin string-like model, the first-order differential operator A(z) writes 

s d A 1 
A{z) = --5-- 


(97) 


(98) 


(99) 


( 100 ) 


dz 2 (z + l) 

The analogue Hamiltonian operators in Eq.s ( |97[ ) and ( |97[ ) form an analogue Supersymmetric 
Quantum Mechanics structure. In fact, defining the analogue charge operators |65j 

0 At 



, Q ] = 


0 0 


( 101 ) 


H = 


( 102 ) 


which are nilpotent, i.e., Q 2 = Q t 2 = 0, and also defining the SUSY-like Hamiltonian 

Hr 0 
0 H l 

we obtain the SUSY-like quantum mechanics algebra |65| 

H = {Q,Q + }, (103) 

and [Q,H] = [Q,H] = 0. The Hamiltonians H R i are related by H l f = Hr and are called 
Hamiltonian superpartners whereas W is known as the superpotential HUES]. 


One remarkable feature of the SUSY-like system (98) is that the massive KK spectrum is the 
same for the both chiralities mmm- Indeed, consider a massive eigenfunction a R (z) of the 


right-handed Hamiltonian H R (97) with mass m R . Define the function 

oil '■= - Aa R . 

m R 


(104) 


Applying the left-handed Hamiltonian H R on the function a^, we find that HlOl{z) = m^d^(z), 
i.e., otL is a left-handed eigenfunction with the same mass of a R . Defining [65] 


a R :=- A^aLi 

m L 


(105) 
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we have Hror^z) = m 2 L aR(z), i.e., olr is a right-handed eigenfunction with the same mass of &l- 
Therefore, for each right-handed eigenfunction d/j exists a left-handed eigenfunction aL with the 
same mass and vice-versa. 

The SUSY-like structure of the Hamiltonians Hrr also guarantees that the spectrum is bounded 
from below. In fact, multiplying any of the Hamiltonians Hrr by the dual eigenfunction d 
respectively, we obtain H-Aor^II 2 = m^joR^H 2 , and hence, m > 0. The absence of tachyonic 
(negative norm) KK modes guarantees the stability of the spectrum. Further, it also enables us to 
employ a probabilistic approach to find the resonant modes, as we will discuss in the next section. 

The Hamiltonian factorization and the absence of negative norm states allow us to reduce the 
problem to find the ground state from a second-order differential equation to a first-order differential 
equation. In fact, for the right-handed massless mode Hrql° r = 0 =>■ H-Ad^H = 0, and thereby, 
Aa R = 0, whereas for the left-handed massless mode oP Ll HlqP l = 0 =>- A^a R = 0. Thus, the 
massless modes a° R L satisfy the equation 

dP R)L ± W (z)a% L = 0 , (106) 


whose solution is given by 


-° = P Tf 0 Z mz')dz' 


a R,L ~ e 


(107) 


By Eq. (107), only one chiral massless mode is normalizable, i.e., localizable. Using the change of 


dependent variable Eq. (90), we obtain the expression (66) for the massles mode. Then, for A > 0, 
only the right-handed massless mode is localized on the brane. 


D. Resonant modes 


In spite of the Kaluza-Klein massive modes are not localized at the brane, some massive states 
can exhibit a relatively large amplitude near the brane [12]. These states, known as resonant 
modes, can be obtained by the quantum mechanical analog structure of the massive modes [TO], 
The resonant modes occur for potentials that exhibit a potential well near the brane and for masses 
m 2 up to the maximum value of the potential barrier mm 


In order to solve the Schroedinger-like equation (95), we need to construct the potential function 


(92). For this, we calculated the warp factors and the gauge angular ansatz in the z -variable from 


the numerical integral of Eq. ([3]) using spline interpolation. We plot in the figures [9] and 10 the 
potential functions for both chiralities. Note that there is a potential well allowing the existence of 
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bound states. The potential has the usual volcano shape. To find solutions of the Schrodinger-like 
equation ( |95| ) with the highest amplitudes near the brane (in comparison with its values far from 
the the defect), we used the resonance method |12j . The relative probability P R .i( rn ) to find a 
particle with mass m in a narrow range 2e around the position 5 of the minimum of the potential 
well may be defined as pans] 


Pr,l ( 7TI ) 


/; max \a R , L (z)\ 2 dz J- z _ 

u min ’ 


rz+e 
J z—e 


\a R:L (z)\ 2 dz, 


(108) 


where z m ; n and z max stand to the domain limits. To perform calculations near the minimum of the 
potential well, we adjusted e = 0.1. 



Figure 9: Potential function for both chirality for 
c = 0.5 and A = 5.0. The thick lines correspond 
to the string-cigar model, while the thin lines, to 
GS one. 



Figure 10: Potential function for both chirality 
for c = 0.5 and A = 9.0. The thick lines corre¬ 
spond to the string-cigar model, while the thin 
lines, to GS one. 


We carried out the numerical integration of the Schrodinger-like equation (95) for the potential 


function (92) into Eq. (108) using the Numerov algorithm [66j for a large sample of the parameters 
c and A. The distribution PrlIjr) exhibited peaks that may be referred as resonant modes 


In the figures 11 and 12, we plot the function Pr,l{™) for c = 0.5. For the left-handed case, there 
are very sharp peaks when A = 4.0 and A = 5.0, while for the right-handed one, when A = 4.0 and 
A = 6.0. However, only the first peak in Pi represents a resonance. To verify this, we solved the 
Schrodinger-like equation for the masses corresponding to each peak in the distribution Pr,l{™)- 


The wave functions are plotted in figures 13 and 14 Note that the solution ai(z) for A = 4.0 has 
the smallest oscillation far from the brane, which characterizes a resonant mode m Although 
the peaks in the distributions -Pfj(m) and P^(m) occur for masses very close (for A = 4.0 when 
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c = 0.5), only the left-handed case has a resonant feature. 

Similar results were obtained for other values of c. In general, the coupling constant A determines 
the existence of a resonant mode, while the geometric parameter c controls the “position” of the 
resonant peaks (i.e. the resonant mass). This is a expected result for a fixed A, since varying the 
geometrical parameter c, which corresponds to change the Planck scale cut-off, different masses 
will be accepted as a resonant state. In a five dimensional thick brane scenario, the resonance of 
fermionic modes was also studied in Ref. [El- 




Figure 11: Plot of the probability distribution Figure 12: Plot of the probability distribution 
PRjJrn) for c: = 0.5 (Left chirality) in the String- Pr,l(to) for c = 0.5 (Right chirality) in the 
Cigar model. String-Cigar model. 


IV. SPIN 3/2 FERMIONIC FIELD 


In this section, we perform the confinement of the Rarita-Schwinger bulk held (spin 3/2). Firstly, 
we start from the following action |32l 3~i1 : 

S& 3/2 = J ^^Mir [M T N T p ^D N ^pd 6 x, (109) 

where square brackets denotes the anti-symmetrization. Its equation of motion has the form 

ri M r N r p ] d n v p = o. (no) 


From now on, we will use the shorter notation y mnp to denote the product of matrices. 

For this spin 3/2 held, the covariant derivative gains an additional term of affine connection 


when compared to spin 1/2 held (56), namely 


Dm^n = ( 9m + ~ i^Am) Hi n ~ ^mn^p- 


( 111 ) 
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Figure 13: Normalized solutions of the 

Schroedinger-like equationin the String-Cigar 
model for the masses corresponding to the peaks 
in the probability distribution. Left and right- 
handed solutions for A = 4.0. 



Figure 14: Normalized solution of the 

Schroedinger-like equation in the String- 
Cigar model for the mass corresponding to the 
second peak in the probability distribution in 
the Fig. ED This is the left-handed solution for 
A = 5.0 and m = 1.387. 


The non-vanishing terms of the Eq. (Ill) with the gauge imposition T# = 


are: 


D^ u = + Jyr (i r r - iqA^j 

IF' ( \F'\ 

D^r = — 'F /i , A-lEV = “ 2 J 

De'k/i = (de + - —T^T r - iqAgj 


( 112 ) 

(113) 

(114) 


Similarly to the decomposition of spin 1/2 in Eq. (59), the Refs. [32., [331 25 J exhibit the 4D 


Rarita-Schwinger vector-spinor in the form 


Vn(x,r,6 ») 



(115) 


Here the spinor 3/2 assumes the KK decomposition: 

1 


^ 4 ) ( x,r,0) = + ^nL nil {x)u Lnl {r) 

v n,l 


aw 


(116) 


where the 4 D section is constrainted by = 0 

m^ v i\) vL , R 


and 7 ^ vp ( d v - iqA v ) ip pR , L = 
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Then, with these restrictions, the non-vanished terms of V MNP Djy'Fp are 


t pixv d^ v = (d„ - iqA tu ) * „ + 


F 

r ppr D^ r = rwDrVp = v™ fa - ^ 

r^A>^ = r p0/i (a* - * g A e ) t /( + ^r^r r ^. 

4xx 


(117) 

(118) 
(119) 


Thus, writing the gamma matrices in the flat form (T 1U = ) and dropping down some 


indexes, for l = 0, the equation of motion (110) with equations (117)-(119) becomes 


d r + 
d r + 


r ) - w ) + W ( r ) U Rni r ) = - 


r ) - w) - W ( r ) U Ln ( r ) = 


URn ( r )> 


( 120 ) 


V F ( r ) 


with T(r) defined in Eq. (64) and W(r) in (|65[). We conclude that the equation (120) is similar to 


the spin 1/2 case presented in the equation 


with the additional term — 


Now, we use the same Ag fixed in Eq. (72), which implies that the W(r) = —AT(r) of Eq. 


(69) holds. Then, in order to obtain a normalized squared modulus solution, the massless-mode of 


rigth-handed spin 3/2 takes the form 

•4(r) = e 0 exp J dr' ((A — 1) 3> — = C 0 


( 121 ) 


where Co is a normalization constant. Afresh, for A = 0, the expression (121) is the same obtained 
in Ref. [32j which is non-normalizable. Moreover, the radial effective action for spin 3/2 is the 
same of the spin 1/2 case presented in Eq. ( |74[ ). Note that there is a correlation between the 
massless modes in Equations (121) and (p7T|) of the form 


l R 


, L (r) = F 2 OL R, L ( r )- 


( 122 ) 


This change promotes a small increase of the amplitude of the zero mode for the spin 3/2. We 
plot in the Figure [15] the massless mode of the Rarita-Schwinger field (121) and compare with Eq. 

©• 

At this point, we have the expressions of the massless modes for the gravitational and scalar 
fields in Eq. (52), 17(1) gauge field in (54), spin 1/2 fields in fflj ) and spin 3/2 in (121). We verified 
that for A = 2 the zero mode of fermionic fields shares similar profile to the bosonic fields. A 


comparative plot is made in Figure 16 
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Figure 15: Comparative plot of right handed 
zero mode of spinorial fields in the String- 
Cigar model. Filled lines corresponds to Rarita- 
Schwinger field, while dashed ones, to spin 1/2 
field. The values for lambda was set to 7 
(broader) and 12 (sharper). The geometric pa¬ 
rameter was set c = 0.5 in both cases. 



Figure 16: Comparative plot of the massless 
modes of bosonic and fermionic fields in the 
String-Cigar model for c = 0.50 and A = 2.0 
(in the fermionic case). The derivative of the 
spinorial fields are indefinite at origin promoting 
similarity with the bosonic fields. 


Using the explicit expressions of the warp factors in Eq. (45) for string-cigar, the zero mode of 
spin 3/2 field can be written as 


u° R (r) = e 0 


(A— 1) 

/ tanh (cr) \ 2 




exp ( -(3 — 5A) [cr + tanh (cr)] 


(123) 


Likewise the spin 1/2 case, for the existing of a zero mode (which vanishes at the origin) is required 
that A > 1, but only for A > 3 there is a normalizable mode with its derivatives null at the origin. 
Moreover, for the left-handed massless modes is required that A —> —A and only one chiral mode 
is allowable to exist in the brane as well. 


A. Spin 3/2 Massive Modes 


Using the conformal radial coordinate z (§, we decouple the equation ( |120[ ) in following two 
second order differential equations: 


F 


aj + 2 :p(.)- ^ I s, + < a =FA)T(z) -a F 




+(i - A 2 )j> 2 o - hLf + ( T 


ur,l = -m 2 u R)L (z), 


(124) 


29 


















where 7(z) = \fFP{z). 


Returning to original variable r, the equation (124) turns to 

"(1^ A) 


dt + 


( 2 / 


+ | ) d r + 


5 f + fg + 10 f + 2 cf + (1 - A 2 ) 


5 f+ 9 

4 J 4 


1 2 


+ 


5 f2 , fg , / 

4 J 4 2 


m 


ur,l(t) = — yU R)L (r). 


( 125 ) 


Minor changes are noted comparing the expressions (125) and (78). In the numerical solution 
of the Sturm-Liouville problem for the Rarita-Schwinger field we verified that the mass spectrum 
is indistinguishable of that presented in Fig. |6j Furthermore, the only difference arises in the 
eingefunctions. The amplitudes for the spin 3/2 are higher than the spin 1/2 case. This can be 
seen in the figure [I7| We intend to perform an analytical study of the spin 3/2 massive modes in 
GS model in a future work. 



Figure 17: Comparative plot of normalized right-handed eigenfunctions for spinorial fields for A = 5.0 and 
c = 0.5 in the String-Cigar model. The most closest mass eigenvalues values was found as mi/ 2 = 0.4024 
and TO 1/2 = 0.4016. 


B. Schrodinger approach for spin 3/2 


Performing a change of the depend variable of Eq. (|124|) in the form 
ur,l(z) = exp 


F 




ur,l(z), 


(126) 


in the Eq. (124), we obtain a Schrodinger-like equation as 


-d't + Vr,£(z) ur,l{z) = m z v,R tL (z), 


(127) 
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where 


V^hz) = A 2 T 2 (z) =f \d z T = W 2 (z) ± d z W. 


(128) 


The spin 3/2 Schrodinger-like equation (128) has the same expression of the spin 1/2 field (92), 
since the change of variables eliminates the multiplicative factor in the first order derivative term 


of ur : l in Eq. (124). As shown in section III C, the form of the potential (128) bears a SUSY-like 


symmetry that guarantees the equality of the KK massive spectra for both chiralities and the 
absence of tachyonic KK states. 


V. CONCLUSIONS AND PERSPECTIVES 

In this work, we have studied the Dirac held in the 6D thick string-cigar braneworld. The string- 
cigar model is a thick string-like scenario which satisfies all the regularity and energy conditions. 
The thin string limit was also analyzed. We considered spin 1/2 and 3/2 (Rarita-Schwinger) bulk 
fields. We made a brief review of the localization of the gravitational and gauge bulk fields. We 
also studied the scalar held and we have shown that it has identical behaviour to that of the 
gravitational held [56] . 

A well-known feature of fermions in a string-like brane is that, unlike the gravitational, scalar 
and gauge helds, it is not possible to trap free fermions with a decreasing warp factor. Due to 
this fact, we proposed a suitable gauge coupling for the fermions with a background gauge held to 
confine both spin 1/2 and spin 3/2 fermions in the brane with positive tension. 

Imposing suitable boundary conditions to guarantee the self-adjointness of the spinor operators, 
a normalized and everywhere well-defined massless mode is obtained for both the thin string and 
string-cigar models depending on the strength of the coupling constant A. This is valid both for 
spin 1/2 as for spin 3/2 case. The massless modes have a shape similar to the energy-momentum 
components of the source of the string-cigar model, whose core is displaced from the origin. This 
shift of the core is a characteristic of string-like branes which source is a vortex for higher winding 
numbers m • Moreover, there is a tiny increase in the amplitude of the Rarita-Schwinger zero 
mode in comparison with spin 1/2 field. 

By means of numerical methods, we have obtained the spectrum and KK eigenfunctions. The 
absence of tachyonic states and the equality of the massive spectrum for the right and the left chi¬ 
ralities is ensured by the supersymmetric Quantum Mechanics structure of the analogue potential. 
For m <C c, the spectrum exhibits the usual behaviour of the Kaluza-Klein theories. Further, we 
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verified that the mass spectrum of the spin 1/2 and 3/2 fields are indistinguishable. For the massive 
modes, the gauge coupling constant plays a similar role with the bulk mass for 5D fermions [63] . 
The massive eigenfunctions behave like those of the thin string-like asymptotically. Near the core 
of the brane, the massive modes are enhanced. This behaviour is also present in the gravitation 
|56j (and, consequently scalar) and vectorial [5?j cases. Further, the amplitudes for the spin 3/2 
are higher than the spin 1/2 case in whole domain. 

Although the KK modes are not localized on the brane, some massive states may exhibit a 
resonant profile. The search of these states was performed by means of the resonance method. 
We found peaks in the probability distribution relating to states which wavefunctions have very 
high amplitudes near the brane. Apart this feature, the oscillation of the wavefunction must be 
as small as possible to characterize a resonant state. The numerical solutions of the Schrodinger- 
like equation showed that this occurs for a specific left-handed mode. Thus, only fermions with 
one chirality is allowed to interact with the brane as a resonant state. An important result is 
that the coupling parameter allows the existence of a resonance peak. The geometric parameter c 
determines the magnitude of the resonant mass which is consistent with the fact that it is related 
to the Planck scale cut-off. Besides, in the Schrodinger approach, the spin 1/2 and 3/2 fields share 
an identical structure, which prevents tachyonic KK modes for both fields. 

Finally, we made a comparative of the normalized massless modes for the bosonic and fermionic 
fields. For a specific value of the gauge coupling, all the held have similar shapes. 

For future works, we intend to study numerically the KK modes, as well as the resonant states, of 
the spin 1/2 and 3/2 fermions in a warped resolved conifold (60| . This scenario provides a geometric 
flow that controls the singularity at the origin in the quantum analogue potential. Moreover, we 
intend to study the properties of gauge coupling term for l / 0 waves and its influence over the 
masseless and the KK modes. 
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